Short-range correlations and shell structure of medium- mass nuclei 
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The single-particle spectrum of the two nuclei 133 Sb and 101 Sn is studied within the framework of 
the time-dependent degenerate linked-diagram perturbation theory starting from a class of onshell- 
equivalent realistic nucleon-nucleon potentials. These potentials are derived from the CD-Bonn 
interaction by using the so-called Vi ow -k approach with various cutoff momenta. The results obtained 
evidence the crucial role of short-range correlations in producing the correct 2sld0g0h shell structure. 

PACS numbers: 21.30.Fe, 21.60.Cs, 21.10.Pc, 27.60.+j 



I. INTRODUCTION 

Shell structure plays a key role in explaining many fea- 
tures of atomic nuclei. In particular, the energy gap be- 
tween groups of levels is crucial in determining their sta- 
bility properties. However, what is the microscopic origin 
of the shell structure is still a subject of great theoretical 
interest. 

In recent years special attention has been focused on 
the role of the various components of the nuclear force in 
the formation and evolution of the shell structure. This 
is mainly related to the advent of new advanced tech- 
niques and facilities, and in particular to the development 
of radioactive ion beams, which allow to study how the 
properties of nuclei may change when going toward the 
proton or neutron drip line. It is for instance a current 
matter of investigation the influence of the tensor force 
on the single-nucleon energies, especially in the case of 
exotic nuclei [E H, S|, where it is expected to produce 
a reduction of the spin-orbit splitting [4]. This aspect 
is also discussed in Ref. Q along with the influence of 
many-body correlations and of states unbound to particle 
emission on the single-particle spectrum of neutron-rich 
nuclei. 

In this work, we have studied the structure of the 
2slG?0g0/i-shell for protons and neutrons, calculating the 
single-particle spectrum of 133 Sb and 101 Sn with realis- 
tic nucleon-nucleon (NN) potentials. Both these nuclei 
are far from stability, 133 Sb having a large neutron ex- 
cess and 101 Sn lying very close to the proton drip line. 
On the other hand, 132 Sn shows strong shell closures for 
both protons and neutrons while this is not the case for 
100 Sn, which may imply different core polarization effects 
on 133 Sb and 101 Sn. We shall see that the softness of the 
core plays an important role in the pcrturbative descrip- 
tion of the single-particle structure of these nuclei. 

The two nuclei are described as a doubly-closed core 
plus one valence nucleon, whose correlations with the 
core are taken into account by means of the time- 
dependent degenerate linked-diagram perturbation the- 
ory Q. In previous works, single-particle valence nuclei 



of the sd-shell have been extensively studied within the 
same framework by other authors 0, H, H|. However, 
the presence of the intruder state Ohu/ 2 makes a more 
challenging problem to reproduce the structure of the 
2slc%0/i-shell. 

We have employed a class of onshell-equivalent realistic 
NN potentials, derived from the CD-Bonn interaction 
[l(| by way of the so-called Vjow-k approach [ll|. This 
approach is a very valuable tool since, by varying the 
cutoff momentum A, it provides different NN potentials 
defined up to A which are all onshell equivalent. This 
allows to study how the short-range (high-momentum) 
components of the NN potential contribute to determine 
the nuclear shell structure. 

The paper is organized as follows. In Sec. II we 
give an outline of our calculations, focusing attention 
on the time-dependent degenerate linked-diagram expan- 
sion. Sec. Ill is devoted to the presentation and discus- 
sion of our results for the 133 Sb and 101 Sn single-particle 
spectra. Some concluding remarks are given in Sec. IV. 



II. OUTLINE OF CALCULATIONS 



As mentioned in the Introduction, we calculate the 
single-particle energies of 133 Sb and 101 Sn within the 
framework of the time-dependent degenerate linked- 
diagram expansion [f| [l2j]. This means that we de- 
scribe the wave function of each nucleus as a single- 
nucleon state whose energy e> is calculated taking 
into account perturbatively the interaction of the odd 
nucleon with the closed core. As a starting point, an 
auxiliary one-body potential U is introduced to break up 
the hamiltonian as the sum of an unperturbed term Hq, 
which describes the independent motion of the nucleons, 
and a residual interaction Hi: 



»=1 »<J 

= (T + (7) + {V NN - U) = H + H x . (1) 



The effective single-particie energies are expressed by 
way of the folded-diagram expansion in terms of the so- 
called Q box: 



e ] =ei+Q-Q' J Q+Q'J Q J Q-Q'J Q J Q J Q+ ... 

(2) 

eg being the unperturbed energy of the single-particle 
state i.e. 



H \j) = (T + U)\j)=el\j) 



(3) 



Eq. ^) represents the well known Kuo-Lee-Ratcliff 
(KLR) folded-diagram expansion [f| [ijj • This expansion 
is written in terms of the vertex function Q-box, which 
is composed of one-body irreducible valence-linked Gold- 
stone diagrams. Once the Q-box has been calculated 
through a certain order in the input interaction Vnn, 
the series of the folded diagrams is summed up to all 
orders using the Lee-Suzuki iteration method 



(a) 



(b) 



FIG. 1: First-order diagram. Graph (a) is the so-called self- 
energy diagram. Graph (b) represents the matrix element of 
the harmonic oscillator potential U = hmtor 2 . 



p , U h h, Hp 





P h 



3* 




FIG. 2: Second-order diagrams. The asterisk indicates non- 
symmetric diagrams, which occur always in pairs giving equal 
contributions. For the sake of simplicity we report only one 
of them. 



Our calculations have been performed including all the 
one-body irreducible valence-linked Goldstonc diagrams 
up to third order in Vnn- I n Fig- Q]thc first-order dia- 
gram, which is composed of the so-called self-energy di- 
agram minus the auxiliary potential [/-insertion, is re- 
ported while the second- and third-order diagrams are 
shown in Figs. and [31 respectively. The [/-insertion di- 
agrams arise in the perturbative expansion owing to the 
presence of the —U term in Hi (see Eq. (JTJ). In our 
calculation we have used as auxiliary potential the har- 
monic oscillator (HO) one, which is the most common 
choice. 

Let us now focus our attention on the input interac- 
tion Vnn- Of course, it would be very desirable to use 
directly a realistic NN potential that reproduces with 
high precision the two-body scattering data and deuteron 
properties. However, to perform nuclear structure calcu- 
lations in a perturbative framework with realistic NN 
potentials one has first to deal with the strong repul- 
sive behavior of such potentials in the high-momentum 
regime. An advantageous method to renormalize the bare 
NN interaction has been proposed in [ll], [3] ■ It consists 
in deriving from Vnn a low-momentum potential Vi ow _k 
defined within a cutoff momentum A. This is a smooth 
potential which preserves exactly the onshell properties 
of the original Vn n and is suitable for being used directly 
in nuclear structure calculations. 

Let us now outline briefly the derivation of Vjow-k [3] . 
The repulsive core contained in Vnn is smoothed by in- 
tegrating out the high-momentum modes of Vnn down 
to A. This integration is carried out with the require- 
ment that the deuteron binding energy and low-energy 
phase shifts of Vnn are preserved by Vi ow _k- This is 
achieved by the following T-matrix equivalence approach. 
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We start from the half-on-shell T matrix for Vnn 
T(k',k,k 2 ) = V NN (k',k) + 



V 



q 2 dqV NN (k',q) 



1 



k 2 — q- 



■T(q,k,k 2 ) , 



where V denotes the principal value and k, kl , and 
q stand for the relative momenta. The effective low- 
momentum T matrix is then defined by 



Ti ow -k(p',P,P 2 ) 



f' x 1 

V \ q 2 dqVi ow - ]l (p',q)— -T low _ k (q,p,p 2 ) 

Jo P -T 



where the intermediate state momentum q is integrated 
from to the momentum space cutoff A and (p 1 ,p) < A. 
The above T matrices are required to satisfy the condi- 
tion 
(4) 

T(p',p,p 2 ) = T low _ k (p',p,p 2 ); (p',p)<A. (6) 

The above equations define the effective low- 
momentum interaction T^ow-ki an d it has been shown 
fllj that their solution is provided by the same KLR 



folded-diagram expansion [6, 12] mentioned before 



(5) 




FIG. 3: Same as Fig. [2] but for third-order diagrams. 
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In addition to the preservation of the half-on-shell T 
matrix, which implies preservation of the phase shifts, 
this Viow-k preserves the deuteron binding energy, since 
eigenvalues are preserved by the KLR effective interac- 
tion. For any value of A, the Viow-k can be calculated 
very accurately using iteration methods. Our calculation 
of Viow-k is performed by employing the iteration method 
proposed in [Til ], which is based on the Lee-Suzuki simi- 
larity transformation [l3l |. 

We have derived from the high-precision CD-Bonn NN 
potential [l(| a class of Vi ow -k's corresponding to val- 
ues of the cutoff A ranging from 2.1 to 2.6 fin -1 . These 
Viow-k's are all onshell equivalent, more precisely they all 
yield the same deuteron binding energy and NN scat- 
tering data up to the anelastic threshold, which corre- 
sponds to a relative momentum k ~ 2.1 fin , as the 
CD-Bonn potential. However, this does not imply that 
these Vjow-k 's give the same results when employed in 
the calculation of the single-particle spectra, since the 
larger the value of the cutoff the more the repulsive high- 
momentum components are explicitly included in Vow-k- 
in other words, Vi ow -k's corresponding to different cut- 
offs have different offshell behaviors (see for istance Cor- 
aggio et al. [IB]). To illustrate this, let us consider the 
offshell tensor force strength. This is related to the in- 
state probability of the deuteron Pd, which implies that, 
when comparing Vow-k's with different cutoffs, offshell 
differences are seen in Pd differences. From inspection of 
Tablefl] where the Pd's for various Vi ow -k's are reported, 
it can be seen that the larger the value of the cutoff the 
stronger the offshell tensor force. 

TABLE I: Values of Pd for Viow-k's corresponding to different 
cutoff momenta A (in fm -1 ) compared with the value yielded 
by the original CD-Bonn potential. 

2.1 2.2 2.3 2.4 2.5 2.6 CD-Bonn 
T96 4~09 4~2l 4~32 4~41 4~49 4^85 



ing to the expression Kuj = 45 A -1 / 3 - 25A -2 / 3 pj} for 
A = 132, and for protons the Coulomb force has been 
explicitly added to the Viow-k potential. 
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FIG. 4: Theoretical single-particle spacings in 133 Sb relative 
to the lds/2 level, as a function of JV max (see text for details). 
The calculations refer to a cutoff A = 2.6 fin. - . The experi- 
mental spectrum of 133 Sb is also reported. 



III. RESULTS 



Calculated and experimental single-particle spacings 
for the one-proton valence nucleus 133 Sb are reported in 
Fig. [5] as a function of the cutoff A. We have found it 
appropriate for the perturbative expansion to limit the 
value of A to 2.6 fin -1 , since larger values of the cutoff 
lead to a rapid deterioration of the the convergence prop- 
erties. The calculations of diagrams have been performed 
using intermediate states whose excitation energy is up 
to iV max = 15 harmonic-oscillator quanta, so to obtain 
the stability the stability of the results within few tens 
of keV when increasing the number of intermediate par- 
ticle states. This is shown in Fig. dj where we report 
the relative SP energies calculated at third order in Hi, 
with a cutoff value A =2.6 fin -1 , as a function of N max . 
The oscillator parameter used is Two = 7.88 MeV, accord- 



It is important to point out that our model space is 
spanned by the 0g 7 / 2 , W5/2, W3/2, 2s 1/2 , and 0/in/ 2 
states above the 132 Sn core, so we have downshifted the 
unperturbed energies of the Ohn/2 and Ogg/2 levels by a 
quantity equal to fku. In this way the 2sld0g0h shell 
is degenerate, as required by the theory of the time- 
dependent degenerate linked-diagram expansion, and the 
O59/2 acts as a hole state below the Fermi surface. Any- 
way, the inclusion of diagrams with (V — £/)-insertions 
in the Q box assures that our results are perturbatively 
independent of the choice of U. In fact, to sum up the 
(V — £7)-insertion diagrams to all orders is equivalent to 
employ a Hartree-Fock self-consistent basis. 
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FIG. 5: (Color online) Calculated and experimental single- 
particle spacings in 133 Sb relative to the ld 5 /2 level (black 
line), as a function of the cutoff momentum A. The blue 
lines are results at third-order in the perturbative expansion. 
The red dashed lines refer to results obtained with the Pade 
approximant [2|1] (see text for details). 



orbit partners have been placed outside the chosen model 
space. The calculated 0g 7 / 2 , 0/i n / 2 relative energies show 
a strong dependence on the cutoff A and a slightly worse 
perturbative behavior. The most striking feature is that 
for small values of the cutoff the discrepancies with re- 
spect to the experimental data are quite large, while a 
better agreement is obtained increasing A. In particu- 
lar, it should be pointed out that with small values of 
the cutoff the 0/in/ 2 level is far away in energy from the 
other levels, thus implying a shell closure at Z = 70. 
A much more reasonable spectrum of 133 Sb is obtained 
when increasing the value of the cutoff which cause the 
intruder Oft.n/2 proton state to join the 50-82 shell. This 
evidences the non-negligible role of the high-momentum 
components of Vnn in the formation of the observed shell 
structure. 

However, after these remarks one may wonder whether 
when increasing A the model space should be en- 
larged to include the 2plf0h levels. More precisely, 
for large value of the cutoff it might occur a sort of 
"shell-quenching" when calculating the single-particle 
spectrum above the closed-shell 132 Sn. As a conse- 
quence, the O/in/2 state as well as its parity partners 
(O/ig/2, 1/7/2, 1/5/2' 2 P3/2, 2p 1/2 ) would approach in 
energy the 2sld0g states, leading to the birth of a sort of 
"macro"-shell 2sld0g2plf0h. 



The results shown in Fig. [5] (blue lines) are normal- 
ized with respect to the 0g?5/ 2 state (straight black line) 
in order to evidence the changes in the spacings when in- 
creasing the cutoff A. In this figure, we have also reported 
(red dashed lines) the results obtained by calculating the 
Pade approximant [2|1] 



[2|1] 



1 - Q3/Q2 



(7) 



Qi being the the vertex function Q-box calculated 
through the ith-order in the inp ut interaction Vnn- Us- 
ing Pade approximants [HI, Il9| | one may obtain an esti- 
mate of the value to which the perturbative series should 
converge. So the comparison between the third-order re- 
sults and those obtained by means of the Pade approxi- 
mant [2|1] is an indicator of the dependence of our results 
on the higher-order perturbative terms. 

From inspection of Fig. [5l it can be seen that the calcu- 
lated spin-orbit splitting between ld 3 / 2 and l<i5/ 2 levels 
is almost independent from A, the experimental sepa- 
ration being reasonably well reproduced. The distance 
between the 2si/ 2 energy and the centroid of the Id lev- 
els is rather independent from the cutoff too, but there 
is no experimental counterpart. It is worth to point out 
that the lrf3/ 2 and 2s!/ 2 relative energies do not exhibit a 
significant dependence on higher-order contributions, the 
difference between the third-order results and those ob- 
tained with the [2|1] Pade approximant being very small. 

The situation is very different for the calculated rela- 
tive energies of the 0<77/ 2 and 0/in/ 2 levels, whose spin- 
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FIG. 6: (Color online) Calculated and experimental single- 
particle spacings in 133 Sb relative to the lds/2 level (black 
line), as a function of the cutoff A. The blue lines denote the 
2sld0g0h levels, while the red one the 2plf0h levels (see text 
for details). 



In order to investigate this question, we have per- 
formed a kind of theoretical "experiment" . We have car- 
ried out calculations where the 2plf0h proton-particle 
unperturbed energies have been considered degenerate 
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with the 2sld0g ones, and then we have turned on the 
interaction H\. The calculated energies relative to the 
0^5/2 level are reported in Fig. [6] It can be seen that 
the spin-orbit splittings of the 2plf0h levels are practi- 
cally independent of the cutoff A, while the 0hu/ 2 level is 
pushed down from its parity partners. In particular, for 
large value of the cutoff the Ohu / 2 state is much closer in 
energy to the 2sld0g levels, confirming the reliability of 
a model space spanned only by the five 2sld0g0h states. 

We have found it interesting to investigate the coun- 
terpart of the 2sld0g0h shell for neutrons, namely the 
single-particle spectrum of 101 Sn. In such a case, the 
closed-shell core is the exotic N — Z nucleus 100 Sn. To 
this end, we have performed the same calculation as for 
133 Sb, except that the harmonic-oscillator parameter is 
hw = 8.55 MeV. 

In Fig. [7] we report the calculated and experimen- 
tal single-neutron energies for 101 Sn relative to the ld 5 / 2 
state as a function of the cutoff A. As can be seen, the 
main difference with respect to the results shown in Fig. 
[S] relies in the convergence properties of the Ogj /2 — ld 5 / 2 
spacing, that is the only spacing observed in 101 Sn. More 
precisely, for the 0177/2 energy level in 101 Sn we have found 
that the third-order contributions of diagrams (27), (32), 
and (33) of Fig. 02 with respect to the second order dia- 
gram (3) of Fig. [21 are larger than what has been found 
for 133 Sb. All these graphs represent the so-called Tamm- 
Dancoff plus Random-Phase Approximation diagrams. 
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FIG. 7: (Color online) Same as in Fig. but for 101 Sn 



The stronger role of higher-order perturbative terms in 
the linked-diagram expansion is probably due to the fact 
that 100 Sn is a much "softer" core than 132 Sn. 

Apart these concernings about the perturbative behav- 
ior of the linked-diagram expansion, our results confirm 
the need of the high-momentum components of the nu- 
clear force in order to improve the description of the neu- 
tron shell structure. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have calculated the single-particle 
spectra of 133 Sb and 101 Sn within the framework of the 
time-dependent degenerate linked-diagram perturbation 
theory. To investigate how the high-momentum com- 
ponents of the NN potential influence the structure of 
the proton/neutron-2sld0g0/i shell, we have employed a 
class of onshell-equivalent realistic nucleon-nucleon po- 
tentials. These potentials have been derived from the 
high-precision CD-Bonn interaction by way of the 
so-called Viow-k approach, varying the cutoff momentum 
from 2.1 up to 2.6 fm _1 . 

The perturbative behavior of the calculated single- 
particle relative energies is quite good for 133 Sb, while 
is less satisfactory for 101 Sn. This may be traced to the 
higher degree of exoticity of 101 Sn, which is very close to 
the proton drip line. 

Our results show that, while the calculated spin-orbit 
splittings are scarcely dependent on the cutoff value, the 
relative energy of the intruder O/111/2 state is strongly 
sensitive to the explicit inclusion of the high-momentum 
components of the NN potential. We have found that 
for small values of the cutoff the O/in/2 level lies far away 
in energy from the other single-particle levels, while the 
correct shell structure is restored when increasing the 
cutoff. 

However, it should be noted that the lack of the high- 
momentum components could be compensated by the in- 
clusion of an effective three-body force [2(1 HH, which, for 
small cutoffs, may play a non-negligible role in determin- 
ing the single-particle spectrum [22j. 
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